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Abstract 

C^ ', By using the density matrix renormalization group (DMRG) technique, the incommensurate 

quantum Frenkel-Kontorova model is investigated numerically. It is found that when the quantum 
Ph ' fluctuation is strong enough, the ^-function featured by a saw-tooth map in the depinned state 
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, ^, . will show a different kind of behavior, similar to a standard map, but with reduced magnitude. 

The related position correlations are studied in details, which leads to a potentially interesting 
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[■"» ■ application to the recently well-explored phase transitions in cold atoms loaded in optical lattices. 
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I. INTRODUCTION 

As a generic model to demonstrate the fractal "devil's staircase" , ID Frenkel Kontorova 
(FK) model has been investigated extensively in the field of classical nonlinear physics 
ever since 1970's.— This model describes the dynamics of a chain of particles connected 
by springs in the presence of a sinusoidal external potential. Its most essential feature is 
the involvement of two competing length scales. One is the average distance between the 
particles and the other is the spatial period of the external potential. This characteristics 
can be demonstrated in various physical systems, such as crystal dislocations,— epitaxial 
monolayers on the crystal surface,— charge density waves^ and dry frictions^ et al. Hence 
studying FK model can lead to the understanding of many universal properties in those 
realistic systems. For example, in the classical regime, when the ratio of two length scales 
are rational, i.e. in the commensurate case, any magnitude of external potential will take the 
particles into pinned state. But for irrational ratio, i.e. in the incommensurate case, only 
when the external potential is above a threshold can the ground state stay in the pinned 
state. This later phenomenon is closed related to the breakup of KAM curves into cantorus 
since the equilibrium positions of the particles are described by the well known Chirikov 
standard map.— This is the so-called Transition of Breaking Analyticity (TBA) pioneered by 
Aubry et al^'^ 

So far, the classical properties of the ID FK model have been well understood. As 
a further step, we need to know how quantum effect will modify those classical results. 
The research along this line will not only bring us into the topic of the classical-quantum 
correspondence, but also take us into a fascinating quantum world, where no classical cor- 
respondence exists. Unfortunately, until now, there is not much work about it because of 
the complexes related to the notorious quantum many-body problems. The first work in 
this direction appeared in 1989^™ in which by using the Feynman path-integral quantiza- 
tion scheme and the Metropolis algorithm, Borgonovi et al. has shown the transition of 
the ^'-function from a standard map to a saw-tooth one for incommensurate quantum FK 
model. Similar results have also been reported by the mean-field theory^ and the squeezed 
state method . ^^i^^i^^ Recently Krajewski et al. has studied the commensurate quantum FK 
model^^ with the algorithm of path-integral molecular dynamics. 

In this paper, instead of making use of the known properties of the classical theory as a 



starting point for the investigations, we attenip to study the quantum FK model directly. 
The employed technique is the density matrix renormalization group (DMRG) method.— 
Since its appearance in 1992,^^ through more than one decade of developments,— DMRG 
has now become the most powerful and efficient numerical method to solve one-dimensional 
quantum many-particle systems. A lot of techniques for the efficiency improvement have 
been put forward. The details of how to apply them upon quantum FK model will be 
presented in Section III. 

Besides the theoretical interest, another motivation for us to carry out this work is the 
rapid experimental development in trapping and controlling Bose-Einstein condensates by 
the optical lattices.— >22i2L2^ Especially, since the seminal experiment by Greiner et al. ob- 
serving the phase transitions in cold atoms from superfiuid (SF) to Mott Insulator (MI),— 
the progresses in this field have been explosive. For the moment, one-, two- and three- 
dimensional quantum systems can be readily created in experiments.— Theoretically, to de- 
scribe the underlying physics of these experiments, the Bose-Hubbard model and its various 
extended versions are widely used .^^'^^'^^ But the intrinsic tight-binding approximations of- 
ten require very strong external potentials. On the other hand, when the external potential 
is very weak, the quantum sine-Gordon model^^ might be an optional model so long as the 
continuum approximation is applicable. Between the two limiting situations, the minimal 
model should be quantum FK model, which is actually a discrete version of the sine-Gordon 
model; hence it goes beyond the limit of continuum or hydrodynamic approximations. It 
is hoped that our work could not only shed some light on the interpretation of the present 
on-going experiments, but also help to provide some new insights in modeling and designing 
the future experiments in wider scope. 

In the following section, the FK model is presented and transformed to the second- 
quantized form suitable for the realization of DMRG algorithm. Section 3 explains in detail 
the DMRG procedures used for quantum FK model. Section 4 is devoted to the results and 
discussions. In Section V, some physics underlying our calculations are further explored. 
Finally we will summarize and conclude the work. 



II. HAMILTONIAN FOR QUANTUM FK MODEL 



Fig. 1 shows the schematic diagram of the system we study. There are A^ + 2 parti- 
cles connected by A^ springs and put in a ID external sinusoidal potential, which could 
be regarded as representing the effect of a loading optical potential. The fixed-boundary 
condition is used with xq = and xn+i = (A^ + 1) a, in which Xi, (i = 0, ■ ■ ■, A^ + 1), is the 
coordinate of the ith particle and a is the average distance between neighbouring particles. 
The Hamiltonian can be expressed as, 
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where go is the wave number of the external potential, 7 the force constant of the springs 
and V the magnitude of the external potential. By introducing the following dimensionless 
parameters. 



Xi = qoXi, Ui = Xi- ifi, H' = —H, /i = qoa, K 
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we can rewrite Eq. (1) as, 
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But Eq. (3) is unwieldly for the realization of DMRG technique. It is necessary to recast it 
into a second-quantized form. For this purpose, the following substitution is employed. 
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in which 6^, 6] are annihilation and creation operators satisfying the boson communication 



relationship, i.e.. 






5ij. Then the Hamiltonian in Eq. (1) is finally reduced to. 



H' = H', + VI, 
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The above Hamiltonian H' looks quite similar in appearance to that of a lattice system 
with bosons tight-bounded to the local potentials and only having the nearest-neighbor 



"hopping" terms. But the physics is quite different. We will come back to this problem 
later in this paper. 

By examining Eq. (5), it can be easily seen that there are three independent tunable 
parameters. The first one is the scaled plank constant h, which is proportional to the 
square root of the particle kinetic energy, E^ ~ fi^qQ/m, in the unit of the particle-particle 
interaction energy 7/^0 > which could be estimated by assuming the particles to occupy the 
lowest energy level in a ID potential with a range around 1/go- ^ embodies the magnitude 
of the quantum influence. Experimentally h can be varied by changing the effective mass m 
of the particle or the effective strength of the particle-particle interactions. 

The second free parameter is K, which measures the strength of the external potential 
in the unit of 7/q'o- This is a parameter most easily accessible in experiments since it is 
directly related to the strength of the loading optical lasers. Generally speaking, the above 
two parameters are of competing effect, namely, h tends to delocalize the particles while 
K to localize them, as is easily understandable and will also be seen from our following 
calculations. 

The last parameter in Eq. (5) is /i = qoa, which denotes the commensurability property 
of the system. Because the commensurate-incommensurate phase transitio n^^i^^ is not the 
focus of this paper. h/2ti is taken to be a fixed irrational number ( -\/5 — 2j /2 in our work, 
as conventionally done in literature. It is already well-known that, in classical regime with 
^ = 0, if yu/27r is a rational number, the system will be in a pinned state for an arbitrary 
external potential. In the incommensurate case, the system can only be in a pinned state 
when K is larger than a critical value Kc] otherwise it will be in a sliding state. But once we 
take quantum effect into considerations, the above scenario will be greatly modified. Before 
going to the details, we will first describe the DMRG procedures used for our calculations 
in the next section. 

III. DMRG PROCEDURE 

The essential steps of DMRG method are to form a superblock consisted of two parts 
denoted by 5* and E. Then, from the targeted state, for example, the ground state |^) of the 
superblock, one can obtain the reduced density matrix of S by tracing out the component 
related to E, i.e., ps = Tte \ip) (^|. Finally by keeping only m eigenstates of ps with the 



largest eigenvalues, all the operators related to S will be rebuilt in the space of reduced 
dimension m. It can be proved that the ground states of the rebuilt S has a maximum 
overlap with the true ground state when S is embedded in the superblock. The error is 
proportional to the sum of the residue eigenvalues thrown away in the rebuilding processes. 
Due to the exponential decrease of the eigenvalues from ps for most of the physical systems, 
very accurate results can still be obtained though working in a much restricted Hilbert space. 
Schematically, the ideas behind DMRG can be explained by Fig. 2. The system S is com- 
posed of left block L and left single site Cl- The environment E is composed of right block 
R and right single site Cr. It is assumed that the Hilbert space of L and R have already 



0f) and |0f),(i = l,...,m). 



been truncated to have dimension m with bases written as _ ^ 
respectively. Initially, the dimension for the full Hilbert space of S is m^. Our aim is 
to seek m most effective bases to truncate this full Hilbert space. DMRG technique has 
provided us an optimal way to realize this aim. For example, we express the targeted 
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state of the superblock as \ip) = Yl Aijki 0f) 0,-^) 
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eigenvectors of the reduced density matrix p^ with largest magnitude of eigenvalues. Here 



^ ^^ are the bases of Cl and Cr. Then the most effective states to represent S is the 
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Ps = J2 ^ijVs^rsfc/ 0i / ^j'^) \4>k^ (W ■ Once the Hilbert space for S is truncated in 

i,j,k,l,r,s=l 

this way, one can add one more single site into S* to make a bigger system and then follow 
the same procedure to make the truncations. Thus step by step, the system size becomes 
bigger and bigger, while the dimension for the effective Hilbert space keeps to be m. This 
is the so-called infinite algorithm in DMRG language. For finite system of fixed number of 
particles, as the size of S grows bigger, E will correspondingly becomes smaller. When S 
reaches the end of the whole system, we will change the role of S and E. The truncation 
procedures are then carried out in the reverse direction along the chain until we get to the 
other end. This kind of forward and backward "sweeping" processes can be repeated many 
times. It is the so-called finite algorithm, which is critical to obtain convergent results in 
DMRG calculations. 

For the Hamiltonian expressed by Eq. (5), because the infinite dimension of the local 
Hibert space, we will not only need to iteratively decimate the space of the growing quantum 
system, but also should make a systematic truncation over the local Hibert space. The 
detailed procedures are summarized as follows. 



1. The Hamiltonian in the local Hilbert space is solved numerically for each single particle 
with the dimension truncated to be M, 
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in which the local basis set for ith particle in the Fock space is expressed as 
\r)^ (r = 1, . . . , M). In coordinate representation, 
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with Hj. to be the Hermite function of order r. Using energy as a weighing mea- 

M 

sure, only states \a)- = J2 Cai |«)j , (a = 1, . . . M), with the lowest m eigen-energy of 
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Eq. (6) are retained to form the working space in the following steps. It should be 
mentioned that when solving H'^^, the following matrix element is used. 
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2. Carry out the finite-size DMRG algorithm with sweeping. The expounding of the 
corresponding technical details in this step can be found in many papers or books.-^^*^ 
We will not repeat it here. Only some specific techniques adopted in our work are 
presented as follows. 

(a) Conventionally, the superblock in DMRG work is composed of two blocks at 
the left and right ends together with two single sites in between, which can be 
schematically expressed as L • •R. If the number of retained states for each block 
is m, the dimension for the superblock would be m^. With m = 20, it would be 



an exceedingly intensive calculation to solve a 160000 x 160000 matrix repeatedly 
in the sweeping process. So we will use L • R scheme as done in Ref. 30 since 
now we only need to solve a 6000 x 6000 matrix for the superblock. The amount 
of calculations are much reduced. 

(b) During the warmup processes with the infinite DMRG algorithm before sweeping, 
the bases for the superblock is selected according to the energy cut-off measure.— 

(c) Once having obtained the solution of the superblock Hamiltonian, the bases of the 
central site are updated by the eigenstates of the density matrix with the neigh- 
boring block parts reduced, i.e., pc = TrL,R\ip) {ipl- If we express the bases of 
left, central and right parts of the superblock by 0f ) , 0f ) , 0f ) , (^ = 1, • • • , '^), 
respectively, according to the decreasing order of the eigenvalues from the cor- 
responding reduced density matrix, we can further reduce the dimension of 



the superblock by keeping only the basis 4>iijk)) ^ w) Vj) 4>k) satisfying 
i + k + j < 171 + 2. Then the superblock matrix size could be decimate again from 
m^ to [{m + 5) (r/i + 4) {m + 3) /6]. For example, if m is taken to be 20, we will 
have a matrix of size 2300 instead of 6000. The effectiveness of this technique 
has been demonstrated in Ref. 32. 

3. After obtaining a convergent result in Step 2, higher energy levels from the local 
Hibert space \a)^ {a = m + 1,M) will be feeded into the systems. To accelerate the 
calculations, we will feed a group of m^ bases each time. Since the added bases might 
not be orthogonal to the old ones for each site, all the bases must be orthogonalized 
and normalized before being submitted to the new round of DMRG sweeping processe 
in the next step. 

4. Continue a new round of DMRG process. But now the central site of the superblock 
L • R will have a dimension of m + iris. As explained in Step 2(b), the same size- 
cut technique for the superblock is used once again. It should be mentioned that 
the procedures adopted in Step 3 and Step 4 is based upon the "optimized phonon 
approach" , which is often used in applying DMRG algorithm for models involving 
phonons.— iS^iS^^ But because Eq. (5) is inhomogeneous along the FK chain, the 
calculations here are much more intensive than those, for example in Holstein model.— 



We have to calculate the effective phonon bases for each particle individually since no 
translation or other symmetries can be used. 

5. Repeat Step 3 by feeding more extra local basis into the systems until we get conver- 
gences in both the decimation of the local Hilbert space and the DMRG sweeping. 

6. Goto Step 1 and repeat all the above procedures to ensure the convergence of the 
results with increasing M. 

In the following calculations, M = 90, m = 10, and rris = 5 are used. The maximum 
dimension for the optimized local Hilbert space is 25, which have been shown to guarantee a 
numerical convergence. For the test running, we have applied our program upon the exactly 
solvable model of a pure harmonic chain.— Very accurate results have been obtained. To 
spare the space, we will not show the details of the test running and the related comparisons 
with the analytical results here. 

IV. RESULTS AND DISCUSSIONS 

A. Hull function and gf-function 

Firstly we check the quantum modifications to the classical hull function and ^'-function, 
which are defined in quantum version^ as, 

Fq{ifi + a) = Xi, 

(y- ^ _ Xi+i-2Xi+Xi_i ^ ^ 

where A denote the quantum average of the operator A over the ground state {ipo)- The 
results are shown in Fig. 3, from which the quantum smoothening effect can be clearly seen. 
Firstly, when h is small, just as expected, there is not much change in the hull and g functions, 
compared with their classical features, i.e., we have the staircases for the hull function and 
the discrete points sitting on a sine-function curve for ^'-function. As h increases until a 
magnitude around 1 ~ 2, the particles will sample more and more classically- forbidden 
points along the chain. At the same time, the ^'-function evolves from a form of standard 
map to a form of saw-tooth map. These results are consistent with those already found in 
literature . ^^i-^^i-^^i^^i^^ From another point of view, they also demonstrates the reliability of 
our numerical codes. 



The most interesting phenomenon takes place in Fig. 3(c) and 3(d), namely, as h increases 
further, the saw-tooth map will become more smoothened and appear in a form similar to a 
standard map but with a reduced magnitude. Actually, this effect has already shown signs 
of its existence in the early work from Borgonovi.— But the fuzzy scattered numerical points 
there deter a clean justification of the obvervance. In our work, because DMRG algorithm 
can delve into very large h, a clear judgement and detailed study are hence possible. In order 
to understand what is the underlying physics for this phenomenon, we recall the squeezing 
state theory,— which tells us, 

9q = ^ «ie 2^"^smXi, (12) 

where Gij represents the correlation function in particle coordinates, i.e., 

no,, = (^ol (x, - Xi) [x, - X,) l^o) . (13) 

As we know, the squeezed state approach meets an almost insurmmoutable challenge in 
calculating Gij directly due to the involvement in seeking the periodic orbits in a 2(A^+ 1)- 
dimension map. But the biggest significance of Eq. (12) lies in the fact that it picks up Gij 
as the most essential physics in understanding the behavior of Qq especially in deep quantum 
regime. Within DMRG framework, Gij can be derived readily. Fig. 4 presents the contour 
map of Gij from which three observances are detailed as follows. 

1. When h is small (<< 1), the correlations are confined to a quite narrow region along 
the diagonal line, which means e~^^'^ ~ 1. So reasonably the system's behavior is 
similar to its classical one according to Eq. (12). 

2. In the region with h around 1 ~ 2, where a saw-tooth map appears in the (^-function, 
we can see some isolated "islands" emerged along the diagonal line in hGij, shown 
in Fig. 4(b). This means some strong quantum correlations and fluctuations among 
the individual parts of the chain. This obvious inhomogeneity of Ga with respect to i 
could be regarded as being responsible for the formation of the saw-tooth map. 

3. As we explore into the deep quantum region with h ^ 2, a regular pattern of Gij can 
be observed, which is quite similar to what happens in a pure harmonic chain. For 
comparisons, we present the case with K = Oin Fig. 4(e). The particles are now highly 
correlated and the previous isolated "islands" merge into each other, which results 
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in a nearly uniform distribution of the diagonal quantum fluctuations. It seems that 
the quantum fluctuations have strongly renormalized the external potential and make 
it almost irrelevant as % flows to infinity. Consequently, according to Eq. (12), the 
homogeneity of Ga will drive the saw-tooth map to the standard-map-like one, but 
with reduced magnitude. 

Now we can produce a physical picture underlying the evolution of Qq with respect Ti. In 
the weak-quantum regime, the particles almost have no correlations. In the middle-quantum 
regime, the particles start to get correlated with each other from some individual parts of the 
systems, characterized by the "islands" in the countour map of Gij. When more and more 
particles becomes correlated, those "islands" will get denser and denser until finally they 
merge into each other and the systems enter a highly- quantum regime. Since the quantum 
correlations are closed related to the tunnelling effect, the specific pattern of Gij in the 
middle quantum regime featured by the saw-tooth map could possibly imply the instanton 
glass phase ever discussed by Zhirov.^° But because our present algorithm is restricted to the 
static ground state, it is very difficult to make the direct investigation of this relationship 
before we can produce the real-time evolution of the systems. Recently the real-time DMRG 
(RTDMRG) simulation has attracted wide attentions and achieved a lot of progresses.— It 
will be our next task to extend the present program to explore the dynamic properties of 
quantum FK model. 

In the above discussions, we have seen that the particle coordinate fluctuation TiGa is a 
very important quantity to characterize the system dynamics under the quantum influences. 
For brevity. Fig. 5 presents directly the variation of UGu with respect to i for different h. 
The "islands" observed in Fig. 4 now becomes peaks. Correspondingly, the special feature 
in the contour map of hGu also appears in Fig. 5. We will not repeat it here. The interesting 
effect needed to be noted in Fig. 5 is that, for % = 15, the curve almost coincide with that 
for a harmonic chain with K = 0. It tells us that in the deep quantum region, the external 
potential plays a negligible role for the system behavior, which corroborates once again the 
reason for the diminishing amplitude in the gf-function. 
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B. Average fixed-distance correlations 

Because there is no translational symmetry in the incommensurate FK systems, the 
following average fixed-distance correlations are defined, 



C. = '-^ . (14) 
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Fig. 6 presents the numerical results of Cm- For comparison, the results for a harmonic 
model without external potential is also plotted. From Fig. 6, we can address three points. 
Firstly, for small h, there are very little correlations between particles, which approaches 
their classical behavior in the limit with h ^ 0. Secondly, with increasing h, longer-range 
particle correlations appear. But instead of monotonous decrease over the particle separation 
as demonstrated by the harmonic model, there appear fiuctuations in the curve. This 
happens just in the parameter range where saw-tooth map effect emerges. Thirdly, as h 
goes beyond about 2, the curves are getting much smoother without any oscillations until 
finally approaches the situation without the external potential. But even in the extremely 
quantum regime, for example h = 15, Cm is still smaller than that in the harmonic model 
except in the short-separation region. This is a bit different from the results shown in Fig. 
5, where nice coincidence of the on-site fluctuations is observed. One reason might be that, 
compared with the accurate calculations of the ground state energy, there is less precision in 
calculating the correlations with DMRG algorithm and more over DMRG normally under- 
estimate the particle correlations. . Another reason is that we are here only considering a 
flnite-size chain. When the distance is as large as the system size, it is understandable to 
expect higher statistical errors for the long-range correlations 

To show the extent of similarity between the FK and the harmonic chains with respect 
to the particle correlations, we introduce the following order parameter, 

N 

L-i ^ m 

C = ^^ , (15) 

E (^Harmonic 
m 
m=l 

which is just the ratio of the sum for the average flxed-distance correlations in FK chain to 
that in the harmonic one. Fig. 7 shows the results, from which we see much clearly that 
as h increases, there is a transition region with h around 1 ~ 2. As /i is large enough, C 
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will saturate. The reason for C to be smaller than 1 should be the same as explained in the 
above paragraph. 

C. Quantum delocalization effect 

Besides the above discussed characteristics of the hull and g'-functions, we have also 
studied the delocalization effect of the quantum modifications based upon the particle wave 
functions. To discuss this properties, the following definitions are used, 

P« (X,) = / ^* (Xi, . . . , X^) V- (Xi, . . . , X;v) n dX,, 

^^^^ " 7}^ P{X)dX ' 

where -ip {Xi, ...,Xj\f) is the many-particle wave functions of the system, P*^*) (Xj) is the 
probabilities to find the ith particle at position Xj and P (X) is the total probability to find 
particles between X and X + dX. Firstly, let us check the localization effect of the external 
potential as done in Ref. 14. Fig. 8 shows the numerical results. It can be easily seen that 
for fixed h, as the external potential becomes deeper, the particles will tend to sit near the 
bottom of the potential until finally all the particles are localized and no particles can be 
found around the top of the potential. This is the localizing effect of external potential. By 
defining an order parameter to describe the probability for particles to sit on the potential 
top, 

Pt = J^llj\i' iX)\'5 [X -2'n[n + - jj dX, (17) 

we can find an abrupt jump in the particle localization as shown in Fig. 9. All the results are 
consistent with those in Hu and Li's paper.— As a further work, we continue to explore into 
the high quantum region by increasing h. Fig. 10 presents the function P{X) by keeping 
K fixed. It is interesting to note that as h increases, more and more particles are found 
away from the potential bottoms. This is understandable since h represents the zero-point 
kinetic energy of the particles. Once again, we plot the disorder parameter Pt in Fig. 11, 
from which we see another kind of abrupt change of P{t) through a transition region with 
h around 1 ~ 2, where the probability to locate the particles on the top of the potential 
suddenly increases. Compared with the localization effect of the external potential, the 
delocalization effect increases much slower. To check the size-dependence of the results, we 
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have shown the calculations with different A^ in Fig. 11, from which no obvious qualitative 
change over the above results can be seen. So the quantum delocalization effects seem 
not much dependent upon the system size. Similar conclusions have also been found in 
calculating other physical quantities. That is why in most of our work, we only confine to 

A^ = 54. 

V. RELEVANCE OF QUANTUM FK MODEL TO OTHER MODELS 

In the above section, we have shown how the quantum effect can greatly modify the clas- 
sical behavior of a FK chain. In this section, we will further discuss the physical significance 
of quantum FK model by exploring its relationship with other models. 

As already mentioned in the introduction, quantum FK model is a minimal model to 
describe a system of cold atoms loaded in an optical lattice between two limiting cases: one 
is the Bose-Hubbard model with tight-binding approximations and the other is the quantum 
sine-Gordon model with hydrodynamic approximations. The marking parameter for these 
models can be represented hj [3 = K/h, which is actually the ratio of the external potential 
to the kinetic energy of the particles. In the following, we will demonstrate how quantum 
FK Hamiltonian Eq. (5) looks like with f3 fiows to +oo and 0. 

Firstly, as /5 ^ +oo, this situation can be realized by taking h ^ with K fixed 
or K —^ +00 with h fixed. Both options mean that the particles are tight-bounded by 
the external potential. Hence the particle index i can be taken as the lattice site index. 
Once this connection is built, we will try to expand H'q in Eq. (5) with respect to h by 
taking ^ ^ 0. To simplify the notation, ^ = \/h/{\^\^) is used for the expansion. By 
utilizing the Baker - Campbell - Hausdorff formula e'^'^^ = e^e^e'i^-^'^^ under the condition 
[A, [A, B]] = [B, [A, B]] = 0, we can have 



cos 



L=o,=oJ'(2i^-j)! 



i {h\ + 6,) + ^^^] = e- V cos (./.) ^ ^ ■ ,,/ ^, , (^1)' (&1 

[ L=Oj=oJ- \^^ - J)- 

oo 2L+1 /I \-f-+l t2L+l 

+-('»EE-|7^^-f-n(^ir(''!) 



The complicated factors in the above sum represents different orders of interactions among 
the local oscillators through the media of the external potential. We will keep only the terms 
up to the order of ^■*. Further more, with random phase approximations (RPA), only the 

14 



operator products fii and n^nj are retained. Then Eq. (5) is reduced to, 






which is stiniulatingly similar in form to the following Bose-Hubbard model with disorders, 

N r TJ -| TV 

H = Y. {ei + ^)h, + -hi{h,-l) -tJ2(brh+i + bl+,k), (20) 

1=1 L ^ -^ i=i 

where t denotes the hopping interaction, U is the on-site repulsion, e^ describes the random 
external potential and /i is the chemical potential. If we further assume no phase coherence 
between b] and bj^i, Eq. (19) and Eq. (20) will appear exactly the same except the fixed 
U in Eq. (20). This implicates that we can use very similar mathematical technique to 
solve both Hamiltonians. But these similarities are only valid literally since the underlying 
physics or interpretations are quite different. Eq. (20) is a true tight-binding single-band 
Hamiltonian. Thus physically, b] |0) represents the creation of a particle on a definite orbit at 
the i-th site. But in Eq. (19), bj |0) denotes driving the i-th. particle onto a local state, \r)-, 
(r = 1, ■ ■ ■, oo), as shown in Eq. (7). In other words, Eq. (19) is more like an infinite-band 
model. Because of this basic difference, we cannot map the results concerning bi, which is 
related to Xi and Pi, directly to the corresponding ones in Eq. (20). 

Now the question becomes, besides as a coordinate operator, what Xi implicates in terms 
of the language from Bose-Hubbard model. The motivation to put forward this question is 
because a lot of work has already done in Bose-Hubbard model. So the answer will surely 
help us to extract more physics from our model by analogy. To get some hint of answering 
this question, we turn to another limit of the quantum FK model with h -^ oo. One obvious 
limit is the pure harmonic chain. In Section V, we have already seen this effect in the study 
of the particle coordinate correlations. 

As is well-known, if only the low-energy properties are considered, a universal effective 
harmonic fluid model can be built for any ID quantum fiuid.^^ For example, if we assume 
the Hamiltonian of a ID bosonic quantum fluid to be, 

r f h^ \ 1 T 

H = dxi/j'^ix) v^ ^(3^) + - / dxdx'V{x - x')n{x)n{x'), (21) 

in which iIj{x) is the annihilation operator at position x satisfying ip{x'),ip''{x) = S{x — x') 
and n{x) = ^"'"(a;)^(a;). Here, it should be noted that, to keep the notation consistent with 
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that in literature, we have dropped the cap over the operator. To avoid confusions, this 
notation is only used in the following discussions of this section. By using 



n{x) = % + n(x). 



(22) 



with uq = N/L the average particle density and assuming that the long-wavelength zero- 
point fluctuations of the density has a wavenumnber -C no, 0(x) and n(x) can be treated 
as conjugate canonical fields. 



[0(x'),n(x)] = -i6{x-x'). 
With these substitutions and approximations, Eq. (21) becomes,— 



H 



£^^^ 



K{d.H^)f + ^{nl[{x'' 



(23) 



(24) 



in which v is the sound velocity and K is a dimensionless parameter determines the quasi- 
long-range order of the bosonic field. When a sinusoidal external potential V{x) = Vcos{qox) 
is added, the above Hamiltonian turns into the famous quantum sine-Gordon form,— 



H 



hv 

2^ 



dx 



K{d,(P{x)f + ^{7rIl{x)f 



-hoo 

+ Vn,JdxcosM{x) + Qx) 



(25) 



where Q is a parameter measuring the commensurability property of the model. By dis- 
cretizing Eq. (25), we obtain a FK Hamiltonian, 



H 



hv 

2^ 



E 



i^[0z-0m]' + ;^(^nO' 



+ Vno ^ cos(go0/ + QO- 



(26) 



Compared with Eq. (1), it is obvious that Pi has mapped into IIj and Xi into 0j. By 
checking Eq. (22) again, we can immediately see that Fig. 4 actually is a display of the 
phase coherence between different particles of a realistic quantum system described Eq. (21) 
in low-energy regime. This is a fascinating scenario since X, and Pi can be readily calculated 
by DMRG algorithm. Since the phase coherence is a purely quantum effect, we can now 
provide another interpretation to the physics in Fig. 4 and Fig. 5. Firstly, when h is 
quite small, the particles almost have no phase coherence. Then as h becomes larger, the 
phase correlation pattern makes the system more like a bose glass, which has been widely 
discussed related to disordered Bose-Hubbard model.i^SdLI^ When h takes the system into 
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highly quantum regime, all particles are in a phase-coherent state, much like a superfluid 
state. With these crude qualitative discussions, we are motivated to follow a remarkable way 
to explore the non-trivial phase diagrams in quantum FK model by investigating the particle 
coordinate and momentum correlations. This will provide us an interesting and attractive 
application of our numerical results to the on-going experiments in the field of cold atoms 
trapped by optical lattices. Further work along this direction will is still in progress and will 
be presented in another paper. 

VI. SUMMARY 

In summary, for an incommensurate quantum FK model in the supercritical regime, we 
can divide the effect of the quantum modifications into three categories. 

1. When h -^1, little alternations can be observed in the hull and g— functions, compared 
with their classical behavior. There is negligible phase coherence among the particles. 

2. As h increases further, for example around 1 ~ 2 when K = 5, individual neighbouring 
particles begin to be highly correlated. During this process, the saw-tooth map will 
appear to represent the ^f- function. 

3. With h high enough, more neighboring particles join to strengthen the correlations 
until all the particles are strongly correlated . At this time, the saw-tooth map will 
be replaced by a map similar to a standard map, but with much smaller amplitude. 

On a qualitative basis, we have shown the relationship of the coordinate correlations 
calculated in our work to the phase coherence properties in a ID bosonic system. The 
interesting variations of the correlation pattern with respect to the quantum modifications 
could have direct relevance to the interference patterns observed in the superfluid/Mott 
phase transitions in cold atoms loaded in optical lattices.— The systematic investigation will 
be the focus of our following work. 
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Figure Captions 

Figure 1 Schematic diagram of FK model with fixed boundary conditions. There are iV+2 
particles connected by springs. The average distance between the neighboring particles is a. 
xq and xn+i are fixed to be and (A^ + l)a. 

Figure 2 Schematic diagram for a superblock composed of left block L, right block R, left 
central site C^ and right central site Cr. L and Cl forms the system part S while R and 
Cr forms the environment part E. 

Figure 3 Variations of the quantum hull function Xi and g functions Qq with respect to 
increasing %. K = 5,N = 56, /i/27r = 34/55. The same set of values for K, N, fi/2ir is used 
for all the other figures except being specially noted. Because the external potential has a 
period 27r, ifi and Xi are folded onto the interval [0, 27r]. 

Figure 4 Contour map of the coordinate correlations HGij between the i-th and the j-th 
particles. For comparison, the case with K = and h = 5 is presented in (e). 

Figure 5 The coordinate fiuctuations hGa of the i-th particle with respect to different 
quantum infiuences under a fixed external potential with K = 5. As in Fig. 4, the results 
corresponding to -ft' = and h = 5 is also presented for comparisons. 

Figure 6 Quantum modifications over the average fixed-distance coordinate correlations 
of the particles Cm under a fixed external potential with K = 5. The case with i^ = ( in 
solid line ) is shown for reference. For clarity, (a) and (b) correspond to different range of 
h. In (a), h < 1 and in (b), h > 1. 

Figure 7 Variations of the order parameter C with respect to h. 

Figure 8 Probability distribution P{X) under different strengths of the external potential 
with h = 0.2. The coordinates have been folded onto the region X G [0, 2tt]. 

Figure 9 Total probability Pt to find the particles on the top of the external potential as 
a function of K for different h. 

Figure 10 Probability distribution P{X) under different strengths of the quantum modi- 
fications. As in Figure 7, the coordinates have been folded onto the region X G [0, 2tt]. 

Figure 11 Total probability Pf to find the particles on the top of the external potential 
as a function of h for different system size A^. 
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FIG. 1: Schematic diagram of FK model with fixed boundary conditions. There are A^ + 2 particles 
connected by springs. The average distance between the neighboring particles is a. xq and xn+i 
are fixed to be and {N + l)a. 
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FIG. 2: Schematic diagram for a superblock composed of left block L, right block R, left central 
site Cl and right central site Cr. L and Cl forms the system part S while R and Cr forms the 
environment part E. 
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FIG. 3: Variations of the quantum hull function Xi and g functions gq with respect to increasing 
?i. K = 5,N = 56, /i/27r = 34/55. The same set of values for K, N, fi/ln is used for all the other 
figures except being specially noted. Because the external potential has a period 2it, i/j, and Xi are 
folded onto the interval [0,27r]. 
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FIG. 4: Contour map of the coordinate correlations hGij between the i-th and the j'-th particles. 
For comparison, the case with K = and ?i = 5 is presented in (e). 
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FIG. 5: The coordinate fluctuations hGa of tlie i-th particle with respect to different quantum 
influences under a fixed external potential with K = 5. As in Fig. 4, the results corresponding to 
K = and ?i = 5 is also presented for comparisons. 
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FIG. 6: Quantum modifications over the average fixed-distance coordinate correlations of the 
particles Cm under a fixed external potential with K = b. The case with i^ = ( in solid line ) is 
shown for reference. For clarity, (a) and (b) correspond to different range of fi. In (a), ?i < 1 and 
in (b), %> I. 
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FIG. 7: Variations of the order parameter C with respect to h. 
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FIG. 8: Probability distribution P{X) under different strengths of the external potential with 
h = 0.2. The coordinates have been folded onto the region X G [0, 2tt]. 
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FIG. 9: Total probability Pt to find the particles on the top of the external potential as a function 
of K for different fi. 
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FIG. 10: Probability distribution P{X) under different strengths of the quantum modifications. 
As in Figure 7, the coordinates have been folded onto the region X G [0, 27r]. 
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FIG. 11: Total probability Pt to find the particles on the top of the external potential as a function 
of ft for different system size A^. 
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